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ABSTRACT 

Many recent observational studies have concluded that planetary systems commonly exist in 
multiple-star systems. At least ~ 20 % of the known extrasolar planetary systems are associated 
with one or more stellar companions. The orbits of stellar binaries hosting planetary systems are 
typically wider than 100 AU and often highly inclined with respect to the planetary orbits. The effect 
of secular perturbations from such an inclined binary orbit on a coupled system of planets, however, is 
little understood theoretically. In this paper we investigate various dynamical classes of double-planet 
systems in binaries through numerical integrations and we provide an analytic framework based on 
secular perturbation theories. Differential nodal precession of the planets is the key property that 
separates two distinct dynamical classes of multiple planets in binaries: (1) dynamically-rigid systems 
in which the orbital planes of planets precess in concert as if they were embedded in a rigid disk, and 
(2) weakly-coupled systems in which the mutual inclination angle between initially coplanar planets 
grows to large values on secular timescales. In the latter case, the quadrupole perturbation from the 
outer planet induces additional Kozai cycles and causes the orbital eccentricity of the inner planet 
to oscillate with large amplitudes. The cyclic angular momentum transfer from a stellar companion 
propagating inward through planets can significantly alter the orbital properties of the inner planet 
on shorter timescales. This perturbation propagation mechanism may offer important constraints on 
the presence of additional planets in known single-planet systems in binaries. 
Subject headings: celestial mechanics - binaries: general - planetary systems 



1. INTRODUCTION 

As of February 2008, among the 220 planet-hosting 
stars (excluding pulsars) 44 (20%) of them are 
memb ers of binary or higher-order multiple-star sys- 
tems (iRashavan et aljj2006|; IDesidera fc: Barbieril [20071 : 
[Eggenberger et al. 2003). Since the systematic adaptive- 
optics searches for stellar companions around known 
planetary systems are not yet complete, the true multi- 
plicit y among planet-hos ting stars should be greater than 
20% fe ggenberger et al.l[2007. ). Many of the stellar com- 
panions around planetary systems are faint, distant, or 
both. Nearly a dozen of them are K-type or later (in- 
cluding eight M dwarf companions and one brown dwarf 
companion). Their sky-projected separations range from 
20 to 12000 AU from the planetary systems, with a me- 
dian of ~ 510 AU. This occurrence of planets in rela- 
tively wide binaries is certainly affected by the observa- 
tional biases. Spectroscopic binaries with angular sep- 
arations 2-6" are often deliberately excluded from the 
planet-search target samples since it is difficult to extract 
precise radial- velocity information fro m the recorded flux 
contr i but ed from the two components ( Valenti fc Fischerl 
[20051 : [Egg enberger et al. 2007). On the other hand, it 
is also shown by theoretical planet formation simula- 
tions that nearby stellar companions may hinder planet 
formation processes in multiple ways. First, a binary 
companion at a distance ~ 20-lOOAU tends to tidally 
truncate the nascent circumstellar disk around the pri- 
mary, th^reb^;_j;edAJ^m^^ 

disk ([Artvmowicz fc Lubow[ [19941 : [Pichardo et~al[ [2005[ : 
[Kiev fc Nelson[[2007[ ). Periodic tidal heating in the disk 
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caused by a companion would also lead to the decay of 
spiral structu r es and evaporation of the volatile materials 
(|Nelson[ [200(1 [Maver et al"][2005[ ). Even if planetesimals 
successfully form in the disk, a close stellar companion 
would stir up their relative velocities so that they are 
likely to fragment through collisions rather than to co- 
agula te ([Kiev fc Nel sonl'2007': 'PaardekooDer,etalJ [200l; 
Theb ault et al.[ [2006: Marzari fc Scholl 20C^^ From 
most of these studies it is inferred that the orbital separa- 
tions of stellar binaries containing planetary systems can 
be preferentially larger than the binary systems without 
planets. 

Distant (and often faint) stellar companions, however, 
may still significantly alter the planetary orbits around 
the pr imary stars on secular timescales. pucker fc Mazeh[ 
(|2002[ ) were the first to investigate the statistical differ- 
ences in the orbital properties of planets in binaries pos- 
sibly caused by such secular perturbations. The most 
distinctive trend they found was the occurrence of mas- 
sive close-in planets exclusively among binaries. Though 
they had only 9 planets in binaries in their samples at 
that time, this trend has been robust for planets with 
msinz > 2Mj and P < 100 days a s the sample size 
increased (see [Udrv fc SantosI [2007[ . and references 
therein). There are also discernible differences in the 
eccentricity-period diagram for different populations of 
extrasolar planets ( Figure [J). It has been pointed out by 
[Eggenberger et all (|2004[ ) that planets in multiple-star 
systems show systematically lower orbital eccentricities 
if their orbital periods are less than 40 days. On the other 
hand, there is a relative paucity of planets in binaries on 
near-circular orbits (e < 0.1), for planets with orbital 
periods greater than 40 days. Indeed, in this region the 
median orbital eccentricity of planets in binaries is 0.36, 
much higher than the median eccentricity 0.27 of planets 



2 



Takeda, Kita & Rasio 



0.6 



0.2 



1^ 



10 AU 
100 AU 
1000 AU 
10000 AU 



'I 



1^ 



• ® ® 

g,® ®# . ^ . 



1 10 100 1000 104 

Orbital Period [days] 

Fig. 1. — Eccentricity-period diagram for 234 extrasolar planets 
(black dots). The red-star symbols represent planets that are as- 
sociated with binary companions, with the symbol sizes scaled to 
1/logp, where p is the projected binary separation. The circled- 
dots are the members of multiple-planet systems. 



around single stars. And lastly, the four largest orbital 
eccentricities and ~ 50 % of the 18 planets with e > 0.6 
are associated with confirmed binary companions. These 
statistical trends strongly suggest that the orbital evolu- 
tion of planetary systems is commonly affected by distant 
stellar companions. 

Theoretical studies have investigated the role of secu- 
lar perturbations by stellar companions in the dynamical 
evolution of planetary systems using simplified models. 
The Kozai mechanism, originally discovered by Kozai 
(fl962) for the orbital analysis of a Jovian asteroid un- 
der the perturbation from the Sun, has been adopted 
for extrasolar planets and successfully modeled some 
of the highly eccentric planetary orbits, assuming that 
those systems co ntain only one planet (|Holman et al.l 
119971 : IWu Mur ray 2003). The Kozai mechanism per- 
mits the planet's orbital eccentricity to stably grow to 
arbitrarily large values up to unity, by tuning the un- 
known mutual inclination between the planetary and bi- 
nary orbits. Recently, there is growing attention on the 
combined effect of the Kozai mechanism and tidal dis- 
sipation mec hanism as an alternative means to produce 
hot J upiters dWu fc MurraA^l2003l : iFabrvckv fc Tremaind 
120071 : IWu et aIT[2007|). For the cases of multiple-planet 



systems in binaries. Marzari e t al.l (1200 5") have systemat- 
ically studied the dynamical outcomes of initially-packed 
triple-planet systems in binaries, for limited cases in 
which the planetary systems are coplanar with the bi- 
nary. 

The secular effect of a non-coplanar stellar compan- 
ion to a coupled system of multiple planets is, however, 
an important problem which has received little attention 
in the past. In reality, planetary systems emerging from 
the nascent circu mstellar disks are more likely to contai n 
multiple planets (|Ida fc Linll2"005l : iThommes et al.ll2007f ). 
It is also rather unlikely that newly f ormed mult iple plan- 
ets orbit in the plane of the binary. iHald (|1994) has sys- 
tematically studied binary systems with solar-type com- 



ponents in search for the alignment of their spin axes, 
by combining the vsini measurements with rotational 
period information. His results showed that spin align- 
ments exist only for relatively close binary systems with 
orbital separations 30-40 AU or less. Because the binary 
separations observed among planet-hosting stars are typ- 
ically much beyond this limit, planetary systems formed 
along the primary's equator are randomly inclined with 
respect to the binary plane. In other words, the dis- 
tribution of the mutual inclination angle / between the 
planetary orbits and the binary plane is uniform in cos/. 
This implies that for about three-quarter of the cases 
planetary orbits are inclined with respect to the binary 
plane by more than 40°, the critical Kozai angle (§ 12. ip . 

The effects of Kozai-type perturbatio ns to multiple- 
planet systems were first investigated by llnnanen et al.l 
fl99 7). While numerically studying the evolution of the 
Solar System under the influence of a hypothetical stel- 
lar companion orbiting at 400 AU from the Sun, they 
discovered that, under certain circumstances, the orbital 
elements of the giant planets evo lve in co ncert, as i f they 
were embedded in a rigid disk. iMalmbe rg et al.J (|2007l ) 
have performed similar numerical experiments but with 
slightly different binary parameters. In their simulations, 
the orbital evolutions of the planets were decoupled, and 
the excited orbital eccentricities of the planets induced 
instability and led to ejection of one or more planets. 

The rigid coupling between planets in resp o nse to the 
binary perturbation found bv llnnanen et al.l (|l997l ) has 
not been understood analytically. A significant challenge 
is involved in developing a generalized perturbation the- 
ory for A/'-body systems with A^ > 3 since the approx- 
imation techniques employed in the traditional secular 
theories are no longer valid. Higher order terms in the 
orbital eccentricities, inclinations, and semi-major axes 
may play important roles in the long-term evolution of 
multiple planets in binaries and thus cannot be elimi- 
nated from the disturbing function as in typical secular 
theories. Systematic mapping of dynamical stabilities us- 
ing numerical integrations is also impractical due to the 
large number of dimensions of the parameter space (the 
mass and the three-dimensional orbits for each planet 
and the binary). 

The objective of this paper is to investigate the secular 
evolution of two planets in wide binaries with a more gen- 
eralized set of initial conditions and to provide analytic 
criteria that separate possible dynamical outcomes. To 
disentangle the complex interplay between binary per- 
turbations and gravitational coupling between planets, 
we first separate the problem into two sets of secular 
three-body problems: the secular evolution of hierar- 
chical, highly inclined triple systems (the Kozai mech- 
anism, § 12. ip and the secular evolution of isolated, near- 
coplanar double-planet systems (the Laplace-Lagrange 
secular theory, § 12. 2p . In § [3l we list the important dy- 
namical classes of double-planet systems in binaries along 
with numerical simulations. The analytic explanation for 
each of the dynamical classes is presented in § [H 

2. SECULAR THREE-BODY PROBLEMS 

In this paper we focus on hierarchical four-body sys- 
tems in which two planets orbit a component of a wide 
stellar binary system. The binary separation is typically 
on the order of 100 AU. The initial planetary orbits are 
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assumed to be nearly coplanar and arbitrarily inclined 
with respect to the binary plane. 

To analyze the long-term evolution of hierarchical four- 
body systems, we first consider the gravitational pertur- 
bation from the stellar companion and the mutual per- 
turbations between the planets separately. The former 
has been studied in detail in the context of the Kozai 
mechanism, which takes account of full three dimensional 
orbits, assuming the hierarchy of the system. The lat- 
ter is the secular evolution of a pair of planets around a 
single star, analytically formulated through the Laplace- 
Lagrange theory, which does not assume any hierarchy 
but uses the assumption of small orbital inclinations and 
eccentricities. Both theories are developed by collecting 
only the secularly varying terms from the time-averaged 
disturbing function. This is an important point, indi- 
cating that these secular perturbations are generally not 
cumulative; a secular perturbation mechanism occurring 
on a longer timescale is averaged out to zero by the 
other perturbation (there are rare cases of secular res- 
onances between two different perturbation mechanisms, 
discussed in § 13. 3p . In the following sections we summa- 
rize the relevant properties of Kozai-type perturbations 
and the Laplace-Lagrange secular theory, and derive the 
secular perturbation timescales of each mechanism. 

Throughout this paper, subscripts 0, 1, 2, and B are 
used to refer to the properties of the primary star, the 
inner planet, the outer planet, and the binary companion, 
respectively. 

2.1. The Kozai Mechanism 

The secular perturbation theory for a hierar chical, 
non-coplanar triple system was first developed by iKozail 
([r962) in the analysis of the motion of a Jovian asteroid 
under the gravitational attraction of the Sun. Here we 
adopt the Kozai mechanism for the orbital evolution of 
a planet around a component of a wide binary system. 

The complete Hamiltonian of the three-body system is 

^_ Gmomi G(mo + mi)mB , 

~ 2ai 2aB • U 

The formulation of the Kozai mechanism takes ad- 
vantage of the hierarchy of the system such that a = 
cii/aB <C 1 and expands the perturbing Hamiltonian 
J^^ in terms of a (since a is small, J^^ is a rapidly con- 
verging series). Th e first terni of J^ ' is the quadrupole- 
order potential JTq (F ord et aD "2QQQV To derive the sec- 
ular equations of motion, short-period terms are elimi- 
nated through averaging the Hamiltonian over the mean 
anomalies , h and Ib- The averaged qu adrupolar Hamil- 
tonian is (jPabrvckv fc Tremain"3l2QQ7[ ). 

- GmomimB aj . . 

mo+mi 8a|(l-e|)3/2 ^ ^ 

X [2 + 3e? - (3 + I2el - lbe\ cos^ uji) sin^ /] .(3) 

Since the binary orbit carries the majority of the angular 
momentum of the system, it is typically taken as the ref- 
erence plane {Ib = 0°) that is invariant in the quadrupole 
approximation. Thus the angle / in Equation [3] corre- 
sponds to the relative inclination of the planetary orbit 
measured from the reference plane. With the averaged 
Hamiltonian, the Lagrange equations of motion can be 



derived (Equations 5 of Innanen et al. 1997) which can be 
numerically integrated for ei(t), /i(t), a;i(t), and Qi{t). 
Now we list the key properties of the Kozai Mechanism. 

(1) Since the mean anomaly li is deliberately removed 
from the averaged Hamiltonian, its canonical conjugate 
momentum Li = mQmi^Gai/ {ttiq + mi) and thus also 
the semi-major axis are constants of motion. Because the 
averaged Hamiltonian is independent of Vti , its canonical 
conjugate i^i, corresponding to the z-component angu- 
lar momentum, is also constant. This implies that the 
averaged Hamiltonian is axisymmetric. Another integral 
of motion, called the Kozai integral, can be defined as 

^ = Jl-el cos/ = Vh. (4) 

The Kozai integral couples the evolution of the eccentric- 
ity and the inclination during Kozai cycles. 

(2) When the mutual inclination between the two or- 
bits is small, the argument of pericenter of a planet sec- 
ularly circulate s in the quad rupolar potential of the bi- 
nary. However. iKozail (|l962l ) discovered that if the inner 
planetary orbit is sufficiently inclined with respect to the 
binary orbit such that h < 0.6, there is a libration solu- 
tion possible for the planetary orbit; the planet's pericen- 
ter argument secularly oscillates around either one of the 
two fixed points, ±90°. If the octupole term is included 
in the perturbing Hamiltonian, the pericenter argument 
in fact is permitted to alternate between circulation and 
libration (Ford et al. 2000). The critical Kozai angle cor- 
responding to h = 0.6 is /koz = 39.23° for a planet with 
initially small orbital eccentricity. 

(3) If the initial mutual inclination Iq is above the crit- 
ical Kozai angle (including a retrograde planetary orbit) 
such that Iq G [/koz, 180° — /koz], then the orbital ec- 
centricity of the planet secularly oscillates with a large 
amplitude. The evolution of the orbital inclination is 
coupled to the eccentricity oscillation through the Kozai 
integral, and / oscillates between /koz and /q. 

(4) The averaged quadrupolar Hamiltonian is con- 
served since it is independent of time. Because the semi- 
major axis is also constant, a conserved quantity can be 
defined from the Equation (|3]), 

C = 2 + 3el — (3 + 12e^ — 15e^ cos^ coi) sin^ / = constant. 

(5) 

This implies that the maximum orbital eccentricity of the 
planet ei^max occurs when coi = ±90°. Thus, ei^max can 
be written as a function of C and /i, which are determined 
from the initial condition. If we assume that the plan- 
etary orbit is initially nearly circular so that ei^o <^ 1, 
then C and h are dependent only on the initial mutual 
inclination /q, and the maximum orbital eccentricity can 
be predicted by a simple formula 



ei,max ^ Y ^ ~ 3 ^^^^ (^) 

This approximation is sufficiently accurate in the limit 
ei,o < 0.2 (Fabrycky & Tremaine 2007). 

(5) The longitude of ascending node of the planet sec- 
ularly circulates, completing one cycle while the orbital 
inclination and the eccentricity oscillate roughly twice. 
This is because of the mirror symmetry of the planetary 
orbit with respect to the / = 90° plane. The rotation 
of the planet's nodal line in the binary plane is equiv- 
alent to and better visualized as the precession of the 
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Fig. 2. — Nodal precession of a planet forced by a binary per- 
turbation. As the nodal line of the planetary orbit rotates in the 
orbital plane of the binary, the angular momentum vector of the 
planet Li sweeps a cone shape around the angular momentum vec- 
tor of the binary Lb • If the initial inclination Iq is above the critical 
Kozai angle /koz? the mutual inclination / also oscillates roughly 
twice between Iq and /koz during one nodal circulation cycle. 

orbital angular momentum vector Li around the normal 
of the reference plane (= the orbital angular momentum 
vector of the binary, Lb). Figure [2] illustrates the nodal 
precession of a planetary orbit in a binary. The angle / 
between the two vectors Li and Lb corresponds to the 
mutual inclination between the orbits. 

The nodal precession of a planet in the Kozai solution 
is always monotonic circulation unlike the pericenter pre- 
cession, which may alternate between libration and cir- 
culation. This is because Qi is measured only from an 
arbitrary reference direction whereas for coi there exists a 
fixed line of reference which is the nodal line of the plan- 
etary orbit in the binary plane. Thus, the circulation of 
the planetary ascending node naturally results from sec- 
ular three-body interaction and it does not require any 
critical angle. 

(6) The characteristic t imescale of the Koza i mecha- 
nism can be estimated as (jKiseleva et al. I [19981 ) 



TKo 



Stt Pi ^ 



2x3/2 mo + mi + mB 

^b) • 

mB 



(7) 



The orbital elements ei, /i, c^i, and fti all evolve roughly 
on this timescale (except that the precession timescale 
of the ascending node is Tkoz)- This is why it is often 
called the "Kozai resonance", althoug h it is not strictly 

Nakailll999r i. 



a resonant interaction (jKinoshita 

(7) When there is another source of secular pertur- 
bation for the planetary orbit on a timescale shorter 
than tkoz, this can suppress the Kozai cycles. Com- 
peting perturbations include general relativistic (GR) 
precession, tidal and rotational bu lges of the host star , 
and other planets in the sy stem (|Wu fc Murravl l2QQ3l : 
iFabrvckv fc Tremaind [20071 ). If any one of these per- 
turbations causes the pericenter of a planet to precess 
faster than does the Kozai mechanism, then the weak 
secular torque from the binary companion secularly av- 
erages to zero, and the Kozai cycles become completely 
suppressed. 

2.2. Secular Gravitational Coupling Between Planets 

The secular evolution of a pair of planets through 
mutual gravitational perturbations has been well- 
formulated in the second-order Laplace-Lagrange secu- 
lar theory ("L-L theory", hereafter). The L-L theory 
also employs orbital averaging and elimination of short- 
period terms to isolate the secular terms in the disturbing 



potential. However, the L-L theory does not assume any 
hierarchy in the planetary system as does the Kozai the- 
ory. The L-L theory is valid for any value of a = ai/a2 
as long as the orbital separation between the planets is 
sufficiently large so that there is no short-term perturba- 
tion. 

In the L-L theory the simplified disturbing function 
for each planet is derived by eliminating the terms that 
are dependent on the mean longitudes (which vary on 
short timescales) and also the terms that are dependent 
only on the semi- major axes (which remain constant). 
Up to the second order in the eccentricities and the in- 
clinations, the eccentricity and inclination perturbations 
are decoupled and thus the disturbing functions can be 
written separately (in the fourth-order terms the evolu- 
tions of_e_and_J_^ecome coupled when they are large; 
see IVeras fc Armit age 2007). It should be also noted 
that the classical L-L theory becomes less accurate in 
the regime with large orbital eccentricities or inclina- 
tions due to the neglected higher-order terms (higher- 
order or semi-analytical secular theories have been re- 
cently formulated to improve the accuracy of the clas- 
sical L-L theory for more general planetary orbits; se e 
iMichtchenko fc Malhotr"all2004l : iLibert fc Henrardll2005[ ) . 
In this paper we employ the L-L theory to derive the ap- 
proximate evolutionary timescales of relevant orbital ele- 
ments of planets. For this purpose, the classical second- 
order L-L theory is still sufficiently accurate, as shown 
in the comparisons with numerical results in the next 
section. 

The second-order L-L theory has its advantage in 
that the secular equations of motion of the planets can 
be easily derived within a simple eigenvalue problem 
(jBrouwer fc Clemencd [19611 : iMurrav fc DermottI fT999h . 
The e-vj part of the disturbing functions yields two 
eigenfrequencies, ^+ and which are functions of 
the planetary masses and semi-major axes (Zhou fc SunI 
2003). The eigenvalues ^+ and g- are the characteristic 
frequencies at which the planets' orbital elements evolve. 
For example, in the secular solutions orbital eccentricities 
of the planets oscillate with a period r = 27r/ — ^_ |. 
The solution for the longitudes of pericenter is more com- 
plex as it is a non-linear combination of the two eigen- 
modes. However, it can be easily shown from the L-L 
solutions that in the following limit 



1 — 3q^/a/b[ 



(1) 

3/2 



< 1 



(8) 



where a = ai/a2, q = m2/mi, and b^^j^^ is the standard 

Lapa lace coefficient (Eq. 7.13 of iMurrav fc DermottI 
Il999f ). the ^_-mode corresponds to the simple pericenter 
precession of the more massive planet, and ^+-mode cor- 
responds to the pericenter precession of the less massive 
planet. Because we are mainly interested in hierarchical 
planetary systems in which the secular perturbation from 
the stellar companion to m2 dominates over the pertur- 
bation from mi to m2, it turns out that ^ is at most 0.1 
or typically much smaller for all the systems discussed 
in § [3 Thus, we will use the derived eigenfrequencies to 
provide estimates for the pericenter precession timescales 
'?"pp,i and Tpp^2 for the inner and outer planets, respec- 
tively, as 
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Fig. 3. — The parameter ^ described in Eq. (O as a function of 
the semi-major axis ratio a for different mass ratios q = vn^lvrix. 
If ^ <C 1, the eigenfrequencies gj^ and g— derived from the L-L 
theory correspond to the pericenter precession frequencies of the 
planets. For all the systems discussed in §[3l ^ <C 0.1. 



PP,i = 27r/^+ 1 
PPJ = 27r/^- J 



with 



(i,j = l,2; i^j). (9) 



Likewise, the 1-Vt part of the disturbing functions 
yields two eigenvalues /i and /2 which describe the char- 
acteristic evolutionary frequencies of the orbital inclina- 
tions and the ascending nodes. Note that there is a subtle 
difference compared to the solution in the e-w space; two 
planets around a spherical star have only one secular fre- 
quency in the space (e.g., /i = 0). This is because 
the orbital inclination is a relative quantity that requires 
an arbitrary reference plane, whereas the orbital eccen- 
tricity is a well-defined quantity and introduces asym- 
metry in the coordinate system. There is in fact only 
one physically meaningful inclination angle for two plan- 
ets, which is their mutual inclination. If, however, the 
oblateness of the star is introduced, it would break the 
degeneracy of the problem, and the L-L theory would 
yield two eigenfrequencies in the I-Q. space. We neglect 
the quadrupolar perturbation from the central star and 
thus derive the single nodal precession timescale rj^pp of 
the planets as 



rnpp = 27r//2 {h 7^ 0). 



(10) 



2.3. Apsidal and Nodal Coupling of Planets 

A pair of planets around a single star can undergo ei- 
ther independent pericenter precessions (apsidal circu- 
lation ) or coupled pericenter precessions (aps idal libra- 
tion) (|Ji et al.l(2003l : iBarnes fc G reenb erdl2006l : iMalhotral 
|2002[ ). If the planets follow apsidal libration, then 
Aw{t) = \zui{t) - zu2{t)\ < Camp + Cfixed at any given 
time t, where Camp < 7r/2 is the libration ampli- 
tude. There are two different modes of apsidal libra- 
tion: Cfixed = corresponds to the aligned libration, and 
Cfixed = TT corresponds to the anti-aligned l ibration. Us- 
ing the secular solutions of the L-L theory, IZhou fc SunI 
(j2003.) have provided the explicit analytical criteria for 



the occurrence of apsidal librations given the initial or- 
bital parameters as 

5 /ia^/' (1 - aV8) 



£2 



initial 



1 — /ia 



1/2 



■cosAtz7o, (11) 



or. 



62 



2 



1 — /ia 



1/2 



1 



initial 5a(l-aV8)cosAti7o 
mi/m2 and Awq = \uj\ 



> 0, (12) 



litial * 



where /i = m^/ m^ ana i\wq = I'^i — '^2|,,,,,,^, 
IZhou fc SunI (|2QQ3f ) called the libration regions defined 
by Equations (pTj) and (p!2|) down- and up-librations, re- 
spectively. 

For two planets orbiting around a single star, whether 
the system secularly undergoes apsidal libration or circu- 
lation can be determined from the initial orbital configu- 
rations. With the perturbation from a stellar companion, 
however, the above criteria need to be evaluated for each 
snapshot of time, and the planets may alternate between 
libration and circulation. 

The role of various external perturbations in af- 
fecting the secular apsidal precession mode of planets 
have been previously investigated by seve ral authors 
(e.g. [Malhotra 2002;' Chiang fc Murray 2002: lFord erall 
l2005f ). Chiang fc Murray (2002) have analyzed the sec- 
ular evolution of the outer two planets c and d of the 
-uAndromedae system and discovered that an adiabatic 
eccentricity excitation to the outermost planet d caused 
by a remnant gaseous disk exterior to the planets would 
naturally bring the planets c and d into small-amplitude 
apsidal libration. In this dynamical scenario, there are 
two effects caused by the outer gas disk to the orbit of 
the outer planet: the secular precession of tu^ (and thus 
change in cos Atuo), and the adiabatic growth of the ec- 
centricity ratio ed/cc- Since 1— /ia^/^ > Ofoi v And c and 
d, these two effects bring the system into the aligned up- 
libration as defined in Eq. (p!2|) . Once this up-libration 
criterion is met, the two planets remain on libration tra- 
jectories, maintaining Aw < 7r/2. As the eccentricity 
growth of the outer planet continues, the libration am- 
plitude further damps down toward Aw = 0. 

Similar apsidal capture mechanism may follow when 
the outer planet of a double-planet system is perturbed 
by a binary companion through the Kozai mechanism. 
Because the condition that the Kozai cycles on the outer 
planet not be suppressed favors orbital configurations 
such that 1 — /iOf^/^ > 0, the most common result from 
a binary perturbation is the capture into aligned up- 
libration (Eq. (p!2|) ). We present specific examples of dif- 
ferent apsidal precession modes as well as their effects on 
the eccentricity evolution of planets in binaries in more 
detail in § 13.31 

The ascending nodes of a pair of planets can also secu- 
larly circulate or librate. The analytical nodal libration 
criteria are (jZhou fc Sunl[2004f ) 



k 

h 



< 



1/2 



initial 



1/2 _ I 



cos Al^o, 



or, 



h 
h 



1 



(13) 



(14) 



initial - COsAfio 

where Af2o = — f^2linitiai- previously mentioned, 
without the presence of external perturbations, the as- 
cending nodes of two planets precess at one characteristic 
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Fig. 4. — The precession of the orbital angular momentum vec- 
tors of double-planet systems in binaries that are (a) in nodal cir- 
culation and (b) in nodal libration. 



frequency. If another source of nodal perturbation, for 
example a secular perturbation from a stellar compan- 
ion, introduces additional nodal precession to one of the 
planets, then the nodal offset Al] of the planets may also 
secularly circulate or librate. 

Figure [4] illustrates the motion of the orbital angular 
momentum vectors for nodally circulating and librating 
systems. In a system in which the nodal coupling be- 
tween the planets is weak compared to an external per- 
turbation, the orbital angular momentum vectors Li and 
L2 independently precess around the normal of the ref- 
erence plane (= Lb for the case of planets in a binary). 
Notice that in a nodally circulating system, the mutual 
inclination angle /12 between the two planetary orbits 
determined as 



- sin /i sin /2 cos (1^1 — 1^2) (15) 



cos /12 = cos 1 1 cos I2 

oscillates within a large range, and the planetary orbits 
no longer remain coplanar. In a nodally librating plan- 
etary systems, the maximum mutual inclination angle is 
limited by the nodal libration amplitude. In Figure[l|(b), 
L2 secularly circulates around Lb due to the perturba- 
tion from the binary while Li precesses around L2 main- 
taining a small angle between the two vectors. If the 
nodal coupling between the planets is sufficiently strong 
and thus the nodal libration amplitude remains small, 
the planetary orbits can remain nearly coplanar as they 
secularly evolve. The rate of nodal precession ~ l/r^^p 
determined from Eq. (p!Q|) is an important measure to 
scale the nodal-coupling strength of planets. If an exter- 
nal perturber precesses the ascending node of one of the 
planets at a greater rate, large-amplitude nodal libration 
or nodal circulation will follow, and the planets will open 
a large mutual inclination angle between their orbits. 

3. DYNAMICAL CLASSES OF DOUBLE-PLANET SYSTEMS 
IN BINARIES 

We find that there are three distinct dynamical classes 
of double-planet systems in binaries: (i) completely de- 
coupled systems in which planetary orbits are indepen- 
dently affected by the binary perturbation; (ii) weakly 
coupled systems in which a large mutual inclination an- 
gle grows between the planetary orbits due to large- 
amplitude nodal libration; (iii) strongly coupled, dynam- 



ically rigid systems in which the inclinations and the as- 
cending nodes of the planets evolve in concert. 

Each dynamical class is presented with a numerical 
simulation. For all the numerical simulation s, we have 
used the integration package MERCURY 6.2 (jChambersI 
I199 9) with the Bulirsch-Stoer mode to accurately ac- 
count for possible occurrence of close encounters between 
planets. The general relativist ic accelerations are also in- 
cluded in the integrations to treat the additional orbital 
precession of planets when they are sufficiently close to 
the primary star. The initial planetary orbits are selected 
to be nearly coplanar (/12 < 5°) and circular (e < 0.1), 
and the Hill stability criterion (|Gladmanlll993[ ) is checked 
to ensure no immediate planet-planet scattering. Each 
simulation is run for at least five Kozai timescales of the 
outer planet (~ 5PI/P2). 

3.1. Evolution of Decoupled Planets in Binaries 

Figure [5] illustrates the orbital evolution of a pair of 
dynamically decoupled planets in a binary following in- 
dependent Kozai cycles. The planetary system is largely 
hierarchical, consisting of an inner Jupiter-mass planet 
at 2AU and an outer ~ IOMq planet at 32 AU. Due 
to the large orbital separation between the planets and 
the small mass of the outer planet, the mutual gravita- 
tional interaction of the planets is completely suppressed 
by the perturbation from the stellar companion. This 
system can be viewed as two non-interacting test par- 
ticles placed in a quadrupole potential of a binary, re- 
sulting in independent Kozai cycles. Since both planets 
are initially inclined by 50° with respect to the binary 
plane, their orbital eccentricities oscillate with roughly 
equal amplitudes. The orbital inclinations of the planets 
Ii and I2 measured from the binary plane also oscillate 
with the same amplitudes, between the initial value of 
50° and the critical Kozai angle /koz ~ 40°. 

Notice, however, that the initial coplanarity between 
the planetary orbits is not maintained at all through- 
out the evolution. Due to the secular perturbation from 
the stellar companion, the orbital angular momentum 
vector L2 of the outer planet precesses on a timescale 
~ 2rkoz,2 ^ 2 Myr, while the orbital angular momentum 
vector Li of the inner planet precesses similarly but on a 
much larger timescale ~ 2rkoz,i ^ 200 Myr (see also Fig- 
ure 2] (a)). Thus the nodal difference /Wt between the 
planets secularly circulates on the timescale rkoz,2, caus- 
ing the mutual inclination I12 (which is a function of Al], 
see Eq. [T5]) to also oscillate on the same timescale with 
large amplitudes. The planetary orbits are maximally in- 
clined with respect to each other every time their nodes 
are anti- aligned, reaching /i2,max ^ 100°, nearly twice 
their initial inclinations with respect to the binary orbit. 

In this example the two planets are dynamically invis- 
ible to each other since their secular mutual interaction 
is completely suppressed by the perturbation from the 
stellar companion. Thus, even though the mutual incli- 
nation I12 between their orbits grows beyond the critical 
Kozai angle /koz for a large fraction of time, there is no 
additional perturbation induced to either planet. 

3.2. Weakly Coupled Planets and Induced Eccentricity 
Oscillations 

Next, we focus on more relevant systems in which (i) 
the evolution of the outer planet m2 is dominated by the 
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Fig. 5. — A pair of gravitationally decoupled planets under- 
going independent Kozai cycles. The initial orbital parameters 
of the system are: for the inner planet (slowly varying black 
curve), mi = l.OMj, ai = 2.0AU, ei = 0.01 and h = 50°; 
for the outer planet (rapidly moving red curve), m2 = 0.032Mj, 
a2 = 31.6 AU, 62 = 0.01 and I2 = 50°; and for the binary com- 
panion, Mb = I.OM0, ttB = 750 AU and cb = 0.20. The two 
ascending nodes are precessing independently. Consequently, the 
mutual inclination I12 (dotted blue curve, scale shown on the right 
axis) between the planetary orbits oscillates on the Kozai timescale 
of the outer planet (rkoz,2 ~ 2 Myr) between 0° and ~ 100°. 




5 10 15 20 

t [Myr] 

Fig. 7. — A weakly coupled, nodally librating double-planet 
system in a binary. The initial orbital parameters of the system 
are: for the inner planet (solid black curve), mi = 0.02Mj, ai = 
0.9 AU, ei = 0.07 and h = 46°, for the outer planet (dashed red 
curve), m2 = 0.6Mj, a2 = 12.7AU, 62 = 0.07 and h = 45°, 
and for the binary companion, Mb = O.7M0, aB = 803 AU and 
cb = 0.78. The binary orbit is initially set at /b =0°. The dotted 
blue curve represents the relative inclination angle /12 between 
the two planets (scale shown on the right axis). The nodal offset 
AQ is librating with large amplitudes. The mutual inclination 
between the planetary orbits grows chaotically, inducing the large- 
amplitude eccentricity oscillations to the inner planet. 




t [Myr] 



Fig. 6. — A weakly coupled, nodally circulating double-planet 
system in a binary. The system is identical to Fig. [5] except that 
the outer planet's mass is increased to m2 = 0.32Mj. Due to the 
increased perturbation by the outer planet, the eccentricity of the 
inner planet oscillates more rapidly compared to Fig. [5] 

secular perturbation from the stellar companion, and (ii) 
the evolution of the inner planet mi is unaffected by the 
stellar companion but affected by the secular torque from 
the outer planet. 

The example in Figure [6] is integrated from the same 
initial conditions as those used in Figure [3 except for 
the outer planet's mass increased by a factor of 10 so 
that m2 = 0.32Mj. The increased perturbation from the 



outer planet now suppresses the binary perturbation on 
the inner planet's orbit. The nodal precession periods 
are ~ 2 Myr for the outer planet (due to the binary) and 
~ 40 Myr for the inner planet (due to the outer planet). 
Even though the evolution of the inner planet is now 
coupled to the outer planet, the planets are in nodal cir- 
culation, i.e., their orbital angular momentum vectors 
Li and L2 precess independently, and the mutual incli- 
nation I12 oscillates between 0° and ~ 100°, similarly 
to the previous example. The difference here, however, 
is that the Kozai cycles of the inner planet are intro- 
duced by the outer planet whose time-averaged inclina- 
tion with respect to the inner planet is ~ 50°. Notice that 
in this case the angular momentum transferred from the 
binary through the outer planet accelerates the eccentric- 
ity oscillation of the inner planet. The Kozai timescale 
of the inner planet is ~ P2 / {Pim2) = 40 Myr, shorter 
than ^100 Myr in the previous example where the inner 
planet is directly perturbed by the stellar companion. 

The secular dynamics of the planets changes drastically 
when the coupling strength between the planets is further 
increased. In Figure [3 the outer gas giant with a mass 
m2 = 0.6Mj at a2 = 12.7 AU follows steady Kozai cycles, 
unaffected by the inner super-Earth with mi = 6AM^ 
at ai = 0.9 AU. The binary perturbation to the inner 
planet is completely suppressed due to the much stronger 
secular torque applied by the outer planet's orbit. 

In this system, the two planets secularly follow large- 
amplitude nodal libration. As illustrated in Figure [41(b), 
in this case the trajectory of the inner planet's orbital 
angular momentum vector Li is a superposition of two 
different motions, namely (i) the precession around the 
binary angular momentum vector Lb , and (ii) the preces- 
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Fig. 8. — A close-in look of the first 7Myr of Figure [3 



sion around the outer planet's angular momentum vector 
L2. The resultant motion of Li is similar to the nutation 
of the spin axis of a gyroscope. When the gravitational 
attraction from the outer planet's orbit is not sufficiently 
strong, Li and L2 can grow a significantly large mutual 
inclination angle /12 (corresponds to the width of the 
cone shape that Li sweeps around L2), as is the case for 
the example in Figure [71 The significant growth of the 
mutual inclination /12 induces additional Kozai cycles to 
the inner planet. Each maximum of ei coincides with 
the minimum of /12, similar to the standard Kozai cy- 
cles (Figure [8]). However, the amplitude of the induced 
eccentricity oscillation is not constant since the angular 
momentum of the perturber m2 is also varying due to 
the interaction with the stellar companion. 

Figure [8] focuses on the first 7Myr of FigureH At t = 
planets are in nearly coplanar, nodally aligned orbits. 
The outer planet's ascending node starts to rotate in the 
plane of the binary due to the perturbation from the 
stellar companion. As the nodal offset AQ between the 
planetary orbits starts to grow, due to the gravitational 
attraction from the outer planet the inner planet begins 
to accelerate its nodal precession. At the same time, the 
orbital inclination /i of the inner planet begins to grow to 
conserve the component of its orbital angular momentum 
parallel to Lb- Shortly before 2Myr, the mutual incli- 
nation I12 between the planetary orbits reaches ~ 40°. 
This initiates the Kozai cycles of the inner planet, and 
ei starts to grow. At t ^ 2.5 Myr the ascending nodes 
of the planets briefiy return to alignment and ei reaches 
the first maximum. At this moment, however, the plane- 
tary orbits are no longer perfectly coplanar because now 
the inner planet's orbit has been further inclined with 
respect to the binary plane. The maximally accelerated 
nodal precession of the inner planet quickly passes ahead 
of the nodal line of the outer planet, thereby growing the 
mutual inclination and decreasing the orbital eccentric- 
ity. When the two nodes are maximally separated around 
t ~ 3.3 Myr, the inner planet's nodal precession changes 
the direction, and ei starts to grow again. Ait ^ 3.6 Myr 
the planetary orbits are again nodally aligned, complet- 




100" 



Fig. 9. — Large-amplitude eccentricity oscillations induced to 
the inner planet of a double-planet system with initial orbital in- 
clination below the Kozai limit /^oz • The initial orbital parameters 
of the system are: for the inner planet (solid black curve) mi = 
0.02Mj, ai = 0.52 AU, ei = 0.09 and h = 14.6°; for the outer 
planet (dashed red curve) m2 = 0.76Mj, a2 = 7.53 AU, 62 = 0.03 
and I2 = 16.5°; and for the binary companion Mb = O.32M0, 
ttB = 419 AU and cb = 0.81. Despite the initially small orbital 
inclination, the eccentricity of the inner planet still oscillates with 
large amplitudes due to the large mutual inclination angle intro- 
duced between the two planetary orbits, while the Kozai cycles on 
the outer planet is significantly suppressed. 



ing one nodal libration cycle. The chaotic eccentricity 
oscillation of the inner planet ensues quasi-periodically, 
corresponding to half the nodal libration period of the 
system estimated from the L-L theory, r^pp/2 ^ 2 Myr. 

We make an additional comment on the effect of the 
GR precession on the evolution of the inner planet seen in 
Figure [71 During the sequence t = 8-14 Myr, ei evolves 
non-periodically, which was not seen in a separate simu- 
lation of the same system without the GR effect. Notice 
that the pericenter argument coi of the inner planet de- 
parts from the ant i- alignment with 002 of the outer planet 
at t ~ 8 Myr as the GR precession accelerates the preces- 
sion of LOi toward apsidal alignment. During this epoch 
of the aligned apsidal libration, the induced Kozai cycles 
on the inner planet are suppressed. After t ^ 15 Myr, the 
system returns to apsidal ant i- alignment, and the chaotic 
oscillations of ei resumes. 

The chaotic eccentricity oscillation of the inner planet 
results from the stellar companion introducing a differ- 
ence in the nodal precession frequencies between the 
planets. As previously mentioned in § 12. 1[ the nodal 
precession of the outer planet naturally follows from the 
secular three-body interactions with the binary stars and 
thus does not require the critical Kozai angle. In the 
sample system shown in Figure [9l the two planets are 
initially inclined by only ~ 15° with respect to the plane 
of the binary, significantly below the critical Kozai an- 
gle. Nonetheless, the nodal precession forced on the outer 
planet's orbit causes the large- amplitude nodal libration 
of the planets, leading to the periodic growth of the mu- 
tual inclination angle. The inner planet's orbital eccen- 
tricity grows to as large as ei ~ 0.8 while the eccentricity 
oscillation of the outer planet is significantly suppressed. 
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Fig. 10. — Evolution of a dynamically rigid double-planet system 
in a binary. The initial orbital parameters of the system are: for 
the inner planet (black solid curve) mi = 0.06Mj, ai = 0.7 AU, 
ei = 0.02 and Ii = 124°; for the outer planet (red dashed curve) 
7712 = 0.22Mj, a2 = 9.1 AU, 62 = 0.09 and h = 129°; and for 
the binary companion Mb = O.93M0, ae = 952 AU and cb = 
0.53. The blue dotted curve shows the relative inclination angle 
Ii2 between the two planets (scale shown on the right axis). In this 
system the ascending nodes of the planets are strongly coupled, and 
the system remains nearly coplanar. The pericenter arguments of 
the planets align briefly near the eccentricity maxima of the outer 
planet, resulting in the synchronous eccentricity oscillations of the 
planets. 



3.3. Dynamically Rigid Systems 

The previous examples highlight the evolution of plan- 
etary systems in which the nodal precession rate of the 
inner planet Qi is small compared to the nodal preces- 
sion rate of the outer planet Ct2 induced by the stellar 
companion. Here we turn our attention to the regime 
Qi > Q2 so that the inner planet can effectively couple 
its nodal evolution to that of the outer planet, thereby 
maintaining a small mutual inclination angle. Figure [TOl 
illustrates the strong nodal coupling of planets. The 
inner Neptune-size planet (mi = 0.06Mj) orbits at 
ai = 0.7 AU, while the outer giant planet (m2 = 0.22Mj) 
orbits at a2 = 9.1 AU. The nodal precession timescale of 
the outer planet is ~ 2rkoz,2 = 20 Myr. The secular nodal 
precession timescale of the inner planet can be estimated 
from the L-L theory and is smaller, r^pp ~ 7 Myr. As 
a result the inner planet successfully follows the nodal 
precession of the outer planet, maintaining tight nodal 
alignment and thus also a small mutual inclination angle. 
Viewed from the orbital plane of the binary, the plane- 
tary orbits precess in concert, as if they were embedded 
in a rigid disk. We call such systems that rigidly respond 
to external perturbations "dynamically rigid" systems. 

As seen in Figure [TOl even though the dynamically- 
rigid planetary orbits remain nearly coplanar, the or- 
bital eccentricity of the inner planet is still periodically 
excited. In this particular system the inner planet's ec- 
centricity oscillation is coupled to the Kozai cycles of the 
outer planet due to the periodic apsidal capture. The 
evolution of the same system in the e2/ei-Aa; space is 
presented in Figure [Til The contours with small dots 
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Fig. 11. — Orbital evolution of the planetary system from Fig- 
ure (TO] presented in the e2/ei-Aa; space. The small dots are the 
secular trajectories, and the blue dashed line is the secular separa- 
trix dividing the regions of libration and circulation, both derived 
from the L-L theory. Different symbols are used to represent the 
different evolutionary stages of the system computed from numer- 
ical integration, with each symbol separated by 10^ yr. The sys- 
tem spends a large fraction of time in the region of near apsidal- 
alignment (small Auj), where the eccentricity ratio 62/61 of the 
planets is nearly fixed. During t ~ 20-40 Myr, the system briefly 
remains in the apsidal libration regime. 



separated with equal time intervals are the secular tra- 
jectories analytically derived from the L-L theory. With- 
out the presence of external perturbations, the planetary 
system would evolve on one of the secular trajectories se- 
lected from the initial conditions. 

The planetary system is initially placed in the secular- 
circulation region, (e2/ei)initiai = 4.5 and Ac^initiai = 
— 117°. As the apsidal offset Au between the planets 
circulates for 360° taking Tpp^i ^ 7 Myr, the system also 
shifts rightward in the diagram to new circulation con- 
tours due to the eccentricity growth of the outer planet 
(and thus growing 62/61) on a timescale rkoz,2 ^ 20 Myr. 
Note that the contours to the right of the phase space are 
progressively more concentrated near the apsidal align- 
ment Alo 0°. This is because the inner planet's peri- 
center precesses much faster when its orbital eccentricity 
is smaller (62/61 greater). Since the eccentricity ratio 
62/61 is roughly fixed to a factor of a few around the 
apsidal alignment, at each episode of apsidal alignment 
the orbital eccentricity of the inner planet is excited. No- 
tice that the eccentricity excitation timescale of the in- 
ner planet is significantly shortened by the outer planet 
bridging the angular momentum exchange from the stel- 
lar companion. Without the presence of the outer planet, 
the eccentricity excitation from the stellar companion 
occurs on a timescale ~ 1 Gyr and thus would be com- 
pletely suppressed by the GR precession. With the outer 
planet aiding the angular momentum exchange between 
the companion star and the inner planet, however, the 
inner planet can evolve on only ~ 20 Myr, the timescale 
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Fig. 12. — A dynamically-rigid planetary system with decoupled 
eccentricity evolutions. The initial orbital parameters of the system 
are: for the inner planet (black solid curve), mi = 0.75Mj, ai = 
1.15 AU, ei =0.04 and h =127°; for the outer planet (red dashed 
curve), m2 = l.SOMj, a2 = 15.1 AU, 62 = 0.10 and h = 129°; 
for the binary companion, Mb = O.8OM0, ae = 1322 AU and 
cb = 0.72. The secular torques are applied to the outer planet by 
the inner planet and also by the stellar companion at similar rates, 
thereby suppressing the amplitude of the outer planet's eccentricity 
oscillation. The inner planet's pericenter argument precesses ~ 10 
times as fast as that of the outer planet, prohibiting the capture 
into apsidal alignment. Consequently, the inner planet's orbital 
eccentricity remains small throughout the evolution. 



of the outer Kozai cycles. 

This periodic apsidal-capture mechanism requires the 
dynamical rigidity of the system and a sufficiently large 
eccentricity oscillation amplitude of the outer planet, as 
well as comparable apsidal precession rates of the plan- 
ets. The sample system shown in Figure [12] does not 
satisfy the last two of these conditions and thus fails to 
excite the inner planet's eccentricity. First, in this ex- 
ample the two competing secular perturbations from the 
inner planet and from the stellar companion both pre- 
cess the outer planet's orbit in ~ 20 Myr. As a result, 
the outer planet's Kozai cycles are suppressed; its eccen- 
tricity reaches only as high as ~ 0.25 while the maxi- 
mum eccentricity predicted from the initial inclination is 
0.58. Also, the apsidal precession timescale for the inner 
planet is r ~ 2 Myr, much shorter than that of the outer 
planet, r ~ 20 Myr, making the episodes of apsidal align- 
ment too brief for the inner planet to adjust its orbital 
eccentricity. 

In an opposite case in which the two apsidal frequencies 
are comparable within a factor of a few, the orbital eccen- 
tricity of the inner planet can be significantly altered at 
each maximum of the outer planet's Kozai cycles. In Fig- 
ure [131 ^1 cind 002 precess due to the outer planet and the 
stellar companion, respectively, both on timescales of a 
few Myr. These similar apsidal precession frequencies al- 
low the long-lasting tight apsidal alignment of the planets 
during each Kozai cycle of the outer planet. Analogous 
to the example in Figure [Til as 62 begins to rise in the 
first 2 Myr, the growing eccentricity ratio 62/61 places the 
system on one of the more elongated circulation contours, 
on which the system quickly returns to the small value of 
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Fig. 13. — Evolution of a dynamically rigid planetary system 
in near secular apsidal resonance. The initial orbital parame- 
ters of the system are: for the inner planet (black solid curve), 
mi = 1.7Mj, ai = 0.7AU, ei = 0.07 and h = 68°; for the outer 
planet (red dashed curve), m2 = 4.5Mj, a2 = 16.4 AU, 62 = 0.03 
and I2 = 65°; and for the binary companion, Mb = O.4M0, 
a-B = 553 AU and cb = 0.37. The pericenter arguments of the 
two planets evolve on similar timescales. At every Kozai cycle 
of the outer planet, the pericenter of the inner planet is strongly 
captured into alignment or anti-alignment with the outer planet, 
resulting in impulsive excitation or damping of ei. 

62/61 thereby steeply exciting 61 at the first Kozai max- 
imum. When the outer planet's eccentricity begins to 
decrease, it occurs on a timescale comparable to the cir- 
culation timescale of cji, thus the inner planet does not 
have enough time to adjust to the circularization of the 
outer planet. Consequently, the eccentricity ratio 62/61 
steadily drops to small values. When the outer planet 
completes its first Kozai cycle, the eccentricity ratio is 
inverted. For the first few Kozai cycles 61 is more effec- 
tively excited than damped, and it systematically drifts 
toward larger values. At the beginning of the fifth Kozai 
cycle of the outer planet {t ^ 11 Myr), however, 62/61 has 
been significantly reduced so that the inner planet gets 
captured into tight apsidal ant i- alignment rather than 
alignment (see Eq. [TT]) . During the ant i- alignment the 
secular torque from the outer planet transfers angular 
momentum to the inner planet's orbit which in turn re- 
duces 61. Consequently, the inner planet's eccentricity 
secularly drifts between and ^ 0.6 as seen in Figure [T3l 

4. DISCUSSION 

4.1. Analytical Boundaries of the Dynamical Classes of 
Planets in Binaries 

Table [U summarizes the coupling properties of the plan- 
ets' orbital elements in each dynamical class. The an- 
alytical boundaries in the parameter space of double- 
planet systems in binaries distinguishing each dynamical 
class presented in the previous section can be provided 
by comparing the several competing secular perturba- 
tion timescales of the systems. The L-L theory yields 
three eigenfrequencies that can be translated into orbital 
evolution timescales Tpp^i, Tpp,2, and r^pp. For a pair 
of planets with a sufficiently large mass or semi-major 
axis ratio such that the outer Kozai cycles are not sup- 
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Fig. 14. — Various secular timescales of double planets in a binary as a function of (a) the orbital semi-major axis and (b) the mass 
of the outer planet, with the following system parameters : mi = 0.3Mj, ai = lAU, uib = IM©, a-B = 500 AU, and cb = 0.5. In the 
left panel (a) m2 = IMj, and in the right panel (b) a2 = 10 AU. In the unshaded area, the Kozai cycles on the outer planet would be 
suppressed by the secular torque from the inner planet. In the light-shaded region, planetary systems undergo dynamically-rigid evolution. 
In the dark-shaded region, the orbital eccentricity of the inner planet grows chaotically due to the growing mutual inclination between the 
planetary orbits. 



TABLE 1 

Dynamical Classes of Double Planets in Binaries 



Dynamical Class 


Description 


AQ 


AI 




ei 


Figure 


Independent Kozai Cycles 


Planets gravitationally decoupled 


Circ. 


Oscillates between and 2/o 


Circ. 


Kozai cycles 


El 


Nodally decoupled 


Accelerated Kozai cycles 


Circ. 


Oscillates between and 2/o 


Circ. 


Kozai cycles ^ 


El 


Large-amplitude nodal libration 


Additional chaotic Kozai cycles ^ 


Lib. 


large 


Lib. ^ 


Kozai cycles ^ 


HE] 


Dynamically rigid 


Synchronized Kozai cycles 


Lib. 


small 


Circ. ^ 


Kozai cycles 


Co] 


Dynamically rigid 


Unaffected ei 


Lib. 


small 


Circ. 


small 


ca 


Dynamically rigid 


Secular ei drift 


Lib. 


small 


Circ. ^ 


large 


ca 



Note. — Circ. and Lib. indicate circulation and libration, respectively. 
^ caused by the outer planet, provided 7o ^ -^koz^ also occurs for Iq < /ko 
caused by the outer planet^ with brief periodic episodes of apsidal capture 

pressed, Tpp^i and Tpp,2 are approximated to the apsidal 
precession timescales of the inner and the outer planets, 
respectively, and r^^pp corresponds to the characteristic 
nodal precession timescale of the planets. A stellar com- 
panion also secularly perturbs the inner and the outer 
planet on timescales rkoz,i and rkoz,2, respectively. Fig- 
ure [HI shows that, generally, for a large region of the 
parameter space Tpp^i < rkoz,i, a regime in which the in- 
ner planet's evolution is affected only by the outer planet. 
For the outer planet, the condition we have imposed is 
such that the perturbation from the stellar companion 
overcomes the perturbation from the inner planet, i.e., 
'?"pp,2 > Tkoz,2, typically corresponding to relatively hier- 
archical systems (the shaded regions in Figure [Mj). 

For such hierarchical systems in which the stellar com- 
panion secularly perturbs the outer planet's orbit, two 
distinctively different evolutionary processes are possi- 
ble, depending on whether the planets can rigidly re- 
spond to the quadrupole perturbations from the stellar 
companion, maintaining a small mutual inclination angle 



^ anti-aligned libration for a large fraction of time with varying amplitude, 

/i2. The key quantity which separates the two dynamical 
classes is the nodal precession rate of each planet. The 
ascending node of the outer planet secularly precesses 
due to the binary companion roughly on a timescale 
2Tkoz,2- The nodal-coupling strength of the inner planet 
to the outer planet can be quantified by the nodal preces- 
sion timescale r^pp derived from the L-L theory. If the 
outer planet's nodal precession occurs faster than that of 
the inner planet such that 

'^Qpp > 2rkoz,2 , (16) 
corresponding to the darker-shaded region of Figure [TH 
then a large mutual inclination angle opens between the 
planetary orbits, and the additional Kozai cycles are in- 
duced to the inner planet on a timescale rind,Koz which is 
shorter than the timescale rkoz,i of a direct perturbation 
from the stellar companion, 

Tind,Koz ~ --p^ < '^koz,2 < ^koz,!- j 

3 Pi 1712 

The amplitude of the induced eccentricity oscillation is 
not constant because the Kozai integral of the inner 
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planet measured in the frame of the outer planet's or- 
bit does not remain constant as the outer planet is also 
subject to angular momentum exchange with the stellar 
companion. 

If, on the other hand, the gravitational nodal coupling 
between the planetary orbits is comparably strong such 
that the nodal precession of the inner planet does not lag 
significantly behind that of the outer planet, i.e., 

TQpp < 2rkoz,2, (18) 

then the planetary orbits can remain close to coplanar. 

In dynamically rigid planetary systems the orbital ec- 
centricity of the inner planet may still evolve, depend- 
ing on the strength of the apsidal coupling between the 
planets. The most extreme case is when the two apsidal 
frequencies are close to the secular resonance such that 

Tpp,i ~ rkoz,2- (19) 

In such cases, the inner planet gets periodically trapped 
into apsidal alignment or ant i- alignment with the outer 
planet. Such repeated apsidal coupling is responsible 
for the secularly-drifting orbital eccentricity of the inner 
planet illustrated in Figure [131 When the pericenter evo- 
lution of the planets is less strongly coupled, the orbital 
eccentricity of the inner planet is periodically excited 
synchronously with the outer Kozai cycles (Figure [TO]). 
If the coupling further weakens such that 

'Tpp,! <C rkoz,2, (20) 

then the two apsidal precessions are completely indepen- 
dent, and the inner planet's eccentricity remains small 
(Figure [E]). 

4.2. Stability of Planets in Binaries 

A stellar companion may induce instability to the 
double-planet system around the primary star that can 
otherwise remain stable for a long term. A possible trig- 
ger of instability is either the Kozai cycles of the outer 
planet, or the chaotic eccentricity oscillation induced to 
the inner planet. In both cases what may follow is the 
loss of a planet either through collision with the host star 
or ejection from the system. The first case of instability 
is rarely observed in our simulations because we have im- 
posed the Hill stability criterion on the initial conditions, 
including the predicted maximum orbital eccentricity of 
the outer planet. The majority of the unstable outcomes 
are found in systems with additional Kozai cycles in- 
duced to the inner planet due to the broken coplanarity 
between the planetary orbits. It needs to be remem- 
bered, however, that our simulations do not include any 
orbital damping mechanism such as tidal dissipation ef- 
fects. Because the eccentricity growth of the inner planet 
induced by the outer planet is rather adiabatic than im- 
pulsive, were tidal effects included in the simulations, 
many of these highly eccentric planets would effectively 
damp their orbital eccentricities and migrate to tighter 
orbits before being lost from the system. 

5. SUMMARY 

We have shown that the complex interplay between the 
secular gravitational perturbations from a stellar com- 
panion and the gravitational coupling between planets 
produces various classes of dynamical outcomes. The two 



distinct evolutionary processes are the dynamically-rigid 
evolution of planetary orbits, and the chaotic growth of 
the mutual inclination angle between the planets. The 
dynamical rigidity of planetary systems can be analyti- 
cally estimated by comparing the nodal precession rates 
of the planets. A weak nodal coupling between the plan- 
ets such that TQpp > rkoz,2 leads to large- amplitude nodal 
libration which steeply inclines the orbital plane of the 
inner planet with respect to the outer planet. Con- 
sequently, the inner planet suffers secular eccentricity 
oscillation with chaotically varying amplitudes. If, on 
the other hand, the gravitational nodal coupling of the 
planets is stronger than the nodal perturbation to the 
outer planet caused by the stellar companion such that 
'TQpp < Tkoz,2, then the orbital planes of the planets pre- 
cess in concert, maintaining near-coplanarity. Even in 
dynamically-rigid planetary systems, the orbital eccen- 
tricity of the inner planet can remain small, or period- 
ically excited in response to the outer Kozai cycles, or 
secularly drift within a large range, depending on the 
apsidal precession frequencies of the planets. 

The periodic eccentricity and inclination perturbations 
from a stellar companion propagating inward through the 
planetary system effectively reduces the orbital evolution 
timescale of the inner planet. Such a cascading angular 
momentum transfer mechanism needs to be incorporated 
to the existing theoretical models of extrasolar plane- 
tary systems in binaries. The traditionally adopted con- 
straint of the critical Kozai angle between the planetary 
systems and the binary orbits may be even relaxed, if 
undetected, non-coplanar planetary companions are ap- 
plying eccentricity perturbations to the detected ones. 
Figure [TH indicates that there is a range of the orbital 
radius of the outer planet within which the additional 
Kozai cycles and significant eccentricity excitation of the 
inner planet are almost certain (a2 ^15 AU for the sys- 
tem presented in Figure [T4|). Also, in a less hierarchical 
system in which dynamical rigidity is ensured (a2 ~ 10- 
15 AU in Figure [T4j), the inner planet's eccentricity can 
be still excited given suitable apsidal-capture conditions. 
Thus, the eccentricity excitation of extrasolar planets by 
binary perturbations is possibly much more efficient than 
previously predicted in the studies b ased on single-planet 
assumptions ( Takeda fc Rasioll2QQ5l ). 

In this paper we have provided a general dynamical 
framework for understanding planetary systems in 
binaries. Yet, there remain some limitations in our 
study which will be further investigated in our future 
papers. First, a broader range of initial planetary 
eccentricities and inclinations needs to be considered. 
Recent theoretical studies have shown that newly- 
formed planets emerging from a protoplanetary disk 
may already possess large orbital eccentricities or 
inclinat ions because of in teract io ns with the gaseous 
disk (|Goldreich fc Saril l2QQ3l : iMoorhead fc Adamsl 
'2007^ or through plan et -planet scattering events 
(Cha tteriee et al.l l2QQ7l : lJuric fc Tremaind l2QQ7l : 
iNagasawa eTall [2008). More realistic distributions 
of initial orbital parameters are necessary to accurately 
determine the likelihood of each dynamical class of 
planets in binaries presented here. Second, we have 
not included tidal dissipation, which may also affect 
the dynamical outcomes of planets in binaries. The 
orbital energy dissipation of highly eccentric planets 
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through tidal friction and subsequent orbital migration 
have been previously stud i ed for hierarchical three-body 
syste n is (IWu fc Murravl l2QQ3l : iFabrvckv fc Tremaind 
I2QQ7I : IWu et all 120071 ). Since significant eccentricity 
growth of the inner planet naturally results when the 
coplanarity between planetary orbits is destroyed, the 
formation of close-in planets among multiple-planet 
systems in binaries through tidal dissipation may be 
fairly common. 

As the observational techniques continue to advance, 
direct measurements of the following observable signa- 
tures will be possible to constrain the dynamical out- 
comes of planets in binaries: (1) the mutual inclina- 
tion angle among multiple-planet systems in binaries 
measured by astrometric observations (e.g., the prelim- 
inary results for ^;Andromedae system by HST Fine 
Guiding Sensor show a large mutual i nclination angle 
~ 35° between the planets c and d; iMc Arthur et al.l 



l2007l ): (2) large spin-orbit misalignments of close-in tran- 
siting planets possibly formed via induced Kozai mi- 
grati on, measured through theRossiter-McLaughlin ef- 
fect (|Fabrvckv fc Tremainl 120071 ): (3) detections of ad- 
ditional planets beyond a few AU from the primary stars 
of binary systems by continuous radial- velocity monitor- 
ing of known planetary systems; (4) detections of stellar 
companions around known hierarchical multiple-planet 
systems by ongoing co-proper motion surveys. 
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